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Abstract 

We study a possible realization of the position- and momentum-correlated atomic pairs that 
are confined to adjacent sites of two mutually shifted optical lattices and are entangled via laser- 
induced dipole-dipole interactions. The Einstein-Podolsky-Rosen (EPR) "paradox" [Phys. Rev. 
47, 777 (1935)] with translational variables is then modified by lattice-diffraction effects. This 
"paradox" can be verified to a high degree of accuracy in this scheme. 

PACS numbers: 03.65.Ud, 34.50.Rk, 34.10.+X, 33.80.-b 



I. INTRODUCTION 

Einstein, Podolsky and Rosen (EPR) Q| put forth the question of whether the quantum 
mechanical description of physical reality is complete, giving the example of a two-particle 
quantum state showing peculiar correlations (dubbed "entanglement" or "Verschrankung" 
by Schrodinger ^): if one measures the position or momentum of one particle, one can 
predict with certainty the outcome of measuring their counterpart for the other particle. 
Thus, depending on which measurement is chosen for the first particle, the value of either 
the position or momentum can be predicted with arbitrary precision for the second particle. 
The ensuing controversy has revolved around the interpretation of the EPR problem and 
its implications on quantum theory ^3']. Later, Bohm considered 0] two entangled spin-1/2 
particles, which have become the focus of attention on this EPR issue: their discrete-variable 
entangled states have served to demonstrate the incompatibility of quantum mechanics with 



local realism, by the violation of Bell's inequality 
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lia ISO, |21|, 122, |23|, 12J, |25|, |26|, 

The original ideal EPR [![ state of two particles — 1 and 2, is, respectively, represented as 
follows in their coordinates or momenta (in one dimension). 



(a;i,X2|V'EPR) = 5{xi - X2), 



(1) 



If two particles are prepared in such a state, and one measures the value of Xi (or pi) 
of particle 1, one can predict the result of measuring X2 (or p2, respectively) with perfect 
precision. The state of Eq. would, however, occupy infinite space and have infinite 
kinetic energy. One can consider more realistic variants of this state, e.g., a Gaussian state 
given by 
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where Ap± = %/Ax± (see Fig. Q). The original EPR state corresponds to the limit 
Ax_/Ax+ 0. After measuring the position of particle 1, the position of particle 2 is 
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FIG. 1: Joint probability distribution of positions (a) and of momenta (6) of EPR-pair ensembles. 
If one measures the position of particle 1, one can predict the position of particle 2 with uncertainty 
~ Aa;_ , whereas if one measures momentum of particle 1, one can predict the momentum of particle 
2 with uncertainty ~ Ap_|_. 
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. In the hmit of Ax_/Ax+ <^ 1, the position of 
particle 2 is centered at Xi, and its uncertainty is ~ Ax_. Similar relations hold also for the 
momentum of particle 2 after the momentum of particle 1 is measured. Thus, either of the 
two conjugate quantities of particle 2 can be predicted with arbitrarily high precision. Of 
course, the Heisenberg uncertainty relation is not violated, since for a single system one can 
measure only one of the two conjugate quantities. 

Approximate versions of the translational EPR state, wherein the (5-function correlations 
are replaced by finite- width distributions, have been shown to characterize the quadratures of 
the two optical- field outputs of parametric downconversion 0,13], or of a fiber interferome- 
ter with Kerr nonlinearity [l9|- Such states allow for various schemes of continuous- variable 
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quantum information processing such as quantum teleportation |2fll . 12 1| or quantum cryp- 



tography 
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large atomic samples 



ar state has also been predicted and reahzed using collective spins of 
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24| . It has been shown that if suitable interaction schemes can be 



realized, continuous-variab 
quantum computation {25 . 



e quantum states of the original EPR type could even serve for 
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Notwithstanding its applications to quantum information processing, the translational 
EPR state of Eq. ((21) does not entail a violation of local realism: such a state has a non- 
negative Wigner function, controlling the position and momentum distribution of each par- 
ticle. Nevertheless, there exist measurement schemes in which an analog of Bell's inequality 
is violated 1^ for such a state — as for any pure entangled state. 

The realization and measurement of the EPR translational correlations of material par- 
ticles appears to be very difficult. There have been suggestions to start with entangled light 
fields and to transfer their quantum state into the state of trapped ions in optical cavities [29| 
or of vibrating mirrors |^ . We have proposed to realize translational EPR states by taking 
advantage of interatom correlations in a dissociating diatom js^. More recently, we have 
considered dipole-dipole coupled cold atoms in an optical lattice as a source of translational 
EPR states 

In order to generate the translational EPR entanglement between interacting material 
particles, one must be able to accomplish several challenging tasks: (a) switch on and off 
the entangling interaction; (b) confine their motion to single dimension, and (c) infer and 
verify the dynamical variables of particle 2 at the time of measurement of particle 1. The 
latter requirement is particularly hard for free particles, sinc e by the time we complete the 
prediction for particle 2, its position will have changed. In we suggested to overcome 
these hurdles by transforming the wavefunction of fiying (ionized) atoms emerging from di- 
atom dissociation by an electrostatic/magnetic lens onto the image plane, where its position 
corresponds to what it was at the time of the diatom dissociation. In [3^ we have proposed 
a solution based on the following steps: (i) controlling the diatom formation and dissociation 
in an optical lattice by switching on and off a laser-induced dipole-dipole interaction; (ii) 
controlling the motion and effective masses of the atoms and the diatom by changing the 
intensities of the lattice fields. In this article we discuss our proposal in more detail and 
elaborate on its principles. 

Our aim here is to demonstrate the feasibility of preparing a momentum- and position- 
entangled state of atom pairs in optical lattices, which would be a variant of the original 
EPR state, owing to lattice diffraction. In Sec. we specify the physical system under 
study. In Sec. Illll the basic properties of single-atom states in optical lattices are discussed. 
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FIG. 2: Proposed scheme of overlapping optical lattices used to prepare the translational EPR 
state. The lattices are displaced from each other in the y direction by They are sparsely 
occupied by two kinds of atoms. Each of the two kinds of atoms feels a different lattice; the shaded 
regions depict the energy minima (potential wells) of the lattices. 



In Sec. II VI we discuss the binding effect of the dipole-dipole interaction. Sec. |V| deals with 
the preparation of EPR states by manipulation of the effective masses of the atoms. In Sec. 
fVTl we discuss experimental demonstration possibilities of measuring the EPR Sec. IVIII is 
devoted to conclusions. 



II. SYSTEM DESCRIPTION 



Let us assume two overlapping optical lattices with the same lattice constant a, as in Fig. 
121 The lattices are very sparsely occupied by two kinds of atoms, each kind interacting with 
only one of the two lattices. This can be realized, e.g., by assuming two different internal (say, 
hyperfine) states of the atoms [33]. In both lattices, the potentials are strongly confining 
in the y and z directions (realized by strong laser fields), whereas in the x direction the 
lattice potential is only moderately to weakly confining. Thus, the motion of each particle 
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FIG. 3: Scheme of the LIDDI interaction: a travehng laser field propagating in the direction x 
along which the atoms are weakly confined. The electric field vector is in the xy plane. The field 
induces dipole moments in the atoms, thereby causing the interatomic interaction. 



is restricted to the x direction. In each direction we assume that only the lowest vibrational 
energy band is occupied. Initially, the potential minima of the lattices are displaced from 
each other by an amount / ^ a in the y direction. An auxiliary laser produces a laser- 
induced dipole-dipole interaction (LIDDI) between the atom pairs. It is linearly polarized 
in the y direction, traveling in the x direction and has a wavelength Ac, moderately detuned 
from an atomic transition that differs from the transition used to trap the atoms in the 
lattice. In the case of two atoms with identical polarizabilities in the geometry of Fig . ^ 
the interatomic LDDI potential induced by a linearly polarized laser is of the form |3J] 



Kid = -VcFeikR), 



(4) 



where 



Fg{kR) = cos {kR cos 9) \ (2-3 cos^ 



cos kR sin kR 
{kR)^ ^ {kR)\ 



+cos^e 
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kR 



and 



(5) 



(6) 



Here the wavenumber is k = 27r/Ac, /c is the coupling laser intensity, and the atomic 
dynamic polarizability a is 



H being the dipole moment element, 00 a the atomic transition frequency, and 00 = kc. The 
position-dependent part Fg{kR) is a function of R, the distance between the atoms, and 
6, the angle between the interatomic axis and the wavevector of the coupling laser. Since 
/ <^ 2a, Vdd(-R) has a pronounced minimum for atoms located at the nearest sites, i? ~ /, 
where 

The LIDDI energy as a function of I and the relative position of the atoms is shown in 
Figs. Eland El Under the above assumptions, we can treat the system as consisting of pairs 
of "tubes", either empty or occupied, that are oriented along x. Only atoms within adjacent 
tubes are appreciably attracted to each other along y, due to the LIDDI. 
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FIG. 4: LIDDI potential as a function of the position of atom 2, given that atom 1 occupies site 
0, for different separations of the two lattices: (a) / = 200 nm, (b) / = 100 nm, (c) / = 40 nm. The 
other parameters are specified in Sec. IVl 
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III. SINGLE- ATOM STATES IN THE WANNIER BASIS 

Let us focus on the subensemble of tube-pairs in which each tube is occupied by exactly 
one atom. In the ID optical lattice, the single-atom Hamiltonian is 

Uq /27rx\ pi 
i^iat = — cos 9 

Here m is the atomic mass, Px is the momentum operator, Uq is the maximum potential 
energy due to the interaction of the atomic dipole with the laser field. 




FIG. 5: Position dependence of the potential energy of atom 2, given that atom 1 is located at site 
0. The horizontal lines denote the lowest band of energies corresponding to uncorrelated atoms 
(bandwidth 4Vhop) and the band of diatom energies Vjd (below the band of uncorrelated atoms). 

where /il is the dipole matrix element of the lattice transition, 5l is the detuning of the 
lattice field from this transition, and II is the intensity of the lattice field. The Hamiltonian 
describes a quantum pendulum. The eigenfunctions of the corresponding Schrodinger 
equation are the Mathieu functions 3^. The eigenvalues form bands, whose spectrum 
depends on the ratio of Uq to the recoil energy. 



Er. 



mXl 



(11) 



so that one can distinguish between strongly binding {Uq ^ -Erec) and weakly binding 
{Uq ~ -Erec) potentials. We assume that the atoms are cooled down to the lowest energy 
band of the lattice, in the absence of LIDDI. 

The state of each atom is then conveniently described in terms of Wannier functions 
\Xj) [36] that are localized at lattice sites labeled by index j. The Wannier functions are 
superpositions of the delocalized Bloch eigenfunctions of the same band, 

1 



\Xj) 



^exp i-ikxj) \(f)k), 



(12) 



where N is the number of lattice sites, and Xj is the position of the jth site. Since the 
Wannier functions are not eigenfunctions of the single-particle Hamiltonian, an atom initially 
prepared in a Wannier state that is localized at one site, will subsequently tunnel to the 
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neighboring sites. Nevertheless, if the tunnehng rate is sufficiently slow, the single-particle 
Hamiltonian eq-Hlat2 in the Wannier basis has a relatively simple form: 
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(13) 



Here the diagonal elements Hq are equal to the energy at the center of the band, and only 
two sets of off-diagonal elements, expressing hopping between the neighboring sites, are 
non- negligible: 



^0 = {Xj\Hle.t\Xj) , Hop = (Xj|Aat|Xj+l), 



(14) 



The hopping rate is related to the e nerg y bandwidth of the lowest lattice band H by 
Vb ~ 4 1 Hop I (for exact expressions see 3q). 

For a moderately deep lattice potential {Uq ^ IS-E'rec), the quantum-pendulum 
Schrodinger equation yields the approximate formulae for Hop and the single-atom effec- 
tive mass: 



Hop 



jEj-ecexp (-0.26-^ 



a?VB 2a2|Vhop| 

The Wannier functions of the lowest band can be approximated by Gaussians: 
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1 3/ i-C -J 



(15) 
(16) 

(17) 



where 
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(18) 



47r2 V 2f/o ■ 

This approximation is relatively accurate for Uo ^ 6-Ercc (see the inset in Fig. fTTjl with 
fidelity > 98%. Note, however, that the Gaussian approximation is not suitable 

for calculating the hopping potential (fT^ since this quantity is very sensitive to the non- 
Gaussian tails of the Wannier wavefunctions. 
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IV. DIATOM BINDING AND TRANSLATIONAL EPR STATES 



Let us now assume that two neighboring tubes in Fig. El are occupied by one atom each 
and the LIDDI is turned on. If the tubes are close to each other (/ <^ Ac, as in Fig. E}, the 
interaction Hamiltonian in the Wannier basis has nonzero elements only for atoms residing 
at the nearest sites, 

^int^KidElxf^)lxf)(xf Kxf I, (19) 
i 

and the total two-atom Hamiltonian is 

This Hamiltonian has been diagonalized numerically, and its eigenvalues are shown in Fig. 
inias a function of the hopping and binding potential strength. One can see that for a suffi- 
ciently large ratio |Vdd|/|Vhop| a band of diatomic states is split off the band of independent 
atoms, towards lower energies. 




IVddl/i^rec IMiopl/^rec 



FIG. 6: Eigenvalues of the two-atom Hamiltonian as a function of (a) the dipole-dipole coupling 
Vdd (for a constant hopping potential |Vhop| = 0.0355 i?rcc), and (6) the hopping potential V|iop (for 
a constant dipole-dipole coupling potential Vdd = 2.16 i^rec)- 

For a strong LIDDI binding, |Vhop| <^ \Vdd\, the ground state of the Hamiltonian eq- 
totalham corresponds to a tightly bound diatom which can be approximated by 

l^o)^^Elxr)lxf)- (21) 
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FIG. 7: Joint probability distribution of the positions of two atoms in the ground state of Hamil- 
tonian ^ for l^hopl = OmBBErec: {a) \Vm\ = l.O^rcc, (6) \Vdd\ = O.lO^rec 

This is a highly correlated state: when particle 1 is found at the jth site of lattice 1, then 
particle 2 is found at the jth site of lattice 2, with position dispersion given by the half- 
width a of the atomic Wannier function in the lowest band, a ^ ctg (Fig. Uk)- The Fourier 
transform of this wave function yields its momentum representation. The corresponding 
momentum probability distribution exhibits anti-correlation similarly to the EPR states 
or (j21), but it reflects the lattice periodicity (Fig. |HK)- In momentum space, the state occupies 
a region of half-width h/{2a), and the probability distribution has narrow ridges along p2 = 
—pi- The width of the ridges is inversely proportional to the lattice size, Ap^ ~ h/{Na), 
and they are shifted by 2TTh/a from each other. 

The probability of atoms to escape their EPR partners "over the next" sites increases 
with the ratio \4iop/Vdd- This leads to an increase of the position dispersion which can be 
estimated by first-order perturbation theory: Let atoms 1 and 2 occupy the jth site in the 
absence of Vhop- With the perturbation \4iop on, atom 2 can occupy also sites j ± 1, which 
have energies |Vdd| above the unperturbed state, with the probability ~ |VhopP/| VddP- This 
contributes to an increase in the diatomic separation dispersion, 



Axi 



Kid/ 



(22) 



resulting in the joint probability distribution of the atomic positions and momenta as a 
function of ^op/Kid) shown in Fig. [7|d and Fig. ISt, respectively. 
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FIG. 8: Joint probability distribution of the momenta of two atoms in the ground state of Hamil- 
tonian ^ for l^hopl = O.OSSS^^rec: {a) \Vaa\ = l-O^rcc, {b) \Vaa\ = O.lO^rec 

The states of the tightly bound diatom form a separate band whose bandwidth is 

Vt'^^^\vilt\ (23) 

below the lowest atomic vibrational band. The diatomic hopping potential V^^^ can be 
estimated by assuming that the two atoms consecutively hop to their neighboring sites, i.e., 
the change 

lxf)lxf)-|xSi)lxfi) (24) 



is realized either via 



xr)lxf)-lxSi)lxf)-lx5i)lxSi), (25) 



or via 



1X^)1X5?!) -lxSr)lxi?r). (26) 
By adiabatic elimination of the higher-energy intermediate states, one obtains 



Vilf - 2^. (27) 

All the states of the diatomic band have correlated positions. However, the momenta are 
not anti-correlated in all these states in the same way as in the diatomic ground state. To 
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realize strong momentum anti-correlations, we have to prepare a state that predominantly 
originates from the bottom of the diatomic band. If we work with thermal states this means 
that the temperature of the system must satisfy 

keT < (28) 

Near the bottom of the band, the diatomic dynamics can be described by means of the 
2-atom effective mass given by 



(2at) 2/i^ ^^iKld 

m — ~ 



(29) 

'a^ ""'hop'-' 

The thermal (kinetic) energy of the diatom is then related to the degree of momentum 
anti-correlation through the sum-momentum spread Ap+ = p^i + Px2, 

Apl ^ ksT2m^^^ ^ ^^^ksT. (30) 

To determine how "strong" the EPR effect is, we compare the product of the half-widths 
of the position and momentum peaks in the tightly bound diatom state with the Heisenberg 



uncertainty limit through the parameter 



31 



' (31) 



2Ax_Aj9+ 

A value of s higher than 1 indicates the occurrence of the EPR effect; the higher the value 
of s, the stronger the effect. 

Strictly speaking, for the multi-peak momentum distribution, one should use a more 
general uncertainty relation, as discussed, e.g., in js?!, that distinguishes the uncertainty 
of multiple narrow peaks from that of a single broad peak. However, even the simple half- 
width of the peaks is a useful measure of the EPR effect. In order to maximize s, we must 
adhere to the trade-off between reducing either Aa;_ , by decreasing |Vhop/Vdd|, or Ap+, by 
increasing |Vhop/Vdd|- The optimum value of s generally depends on the temperature of the 
diatom, as detailed below. 



V. EPR STATE PREPARATION 



Cooling down the diatomic system to prepare the EPR state is a non-trivial task. We sug- 
gest a "cooling" procedure which takes advantage of the difference between the single-atom 
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and the diatom bandwidths, and of the possibihty to change the hght-induced potentials. 
The key is first to cool down individual atoms and then separate the unpaired atoms from 
the diatoms. The scheme consists of three steps: 

(i) We first switch on only an external, shallow, harmonic potential in the x direction (all 
other potentials being off), and cool the x-motion of the atoms down to its ground state. 
The width of the ground state should be several times the lattice constant; it is related 
to the desired momentum anti-correlation by aE ~ h/ {\^Ap^). The temperature necessary 
to achieve this must be 

T ^h^/{AmkBal). (32) 

(ii) A weak lattice potential in the x-direction is then slowly switched on, so that the 
state becomes 

(E«.lxr))(E«dxP^)) = i:o?MT)\xf)+i:<^Mxf)\x?), (33) 
where the coefficients 

a, ~exp[-(j-jo)V/(4a|)] (34) 

are Gaussians localized around the minimum of the external potential. 

(iii) We switch on the LIDDI and change the external potential, from an attractive well to 
a repulsive linear potential, acting to remove the particles from the lattice. The two parts of 
the wavefunction 33 will respond differently: The motion of the paired atoms (corresponding 
to J2j '^'j\xf^)\xf^)) will remain in the vicinity of the initial position because of their narrow 
energy band (Fig. Single (unpaired) atoms, whose bandwidth is substantially larger, 
will travel a much longer distance before hitting the top of the energy band. Thus, after a 
properly chosen time, the unpaired atoms could be removed from the x-region of interest, 

(2,') _____ 

which is ~ Vhop/V^op longer for single atoms than for diatoms [see Eq. (j2Sl)]- 

Provided the time tc is short enough for the system to remain near the bottom of these 
two bands, the dynamics can be interpreted in terms of the appropriate effective masses: 
the "heavy" diatoms with mass move much slower than the "light" single atoms with mass 
(I16|l . Thus, after changing the sign of the external potential, the unpaired atoms will be 
ejected out of the lattice and separated from the diatoms as glumes from grains. 
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FIG. 9: Separating single (unpaired) atoms from diatoms: an external repulsive potential causes 
both the single atoms and the diatoms to move on a surface of constant energy. The diatoms hit 
the top of the energy band after being displaced by a much shorter length than the unpaired atoms. 

This effect is illustrated by the numerical simulation in Fig. for two lithium atoms in 
two lattices with Al = 323 nm (corresponding to the transition 2s-3p) and a dipole-dipole 
coupling field of Ac = 670.8 nm (transition 2s-2p). The dipole moment element of the lattice 
transition is 1.26 x 10^^^ Cm, while the LIDDI coupling dipole element is 2.7 x 10^^^ Cm. 
From these values we get the recoil energy i^rec = 1-85 x lO^^^J. The lattice and LIDDI 
field intensities are Jl = 0.186 W/cm^ and Iq = 0.023 W/cm^. The corresponding field 
detunings are 5l = 507l, 5c = IOO7C, the respective decay rates being 7l = 1.2 x 10^ 
s~^, and 7c = 3.7 x 10^ s~^. The two lattices are displaced by / = 40 nm. From these 
values we get the lattice potential Uq = 3.93-E'rec, the LIDDI potential of the nearest atoms 
Vdd = — 0.5-Erec, and the hopping potential Vhop = — 0.09i?rec- The two-particle hopping 
potential is then V^^^^^ ~ — 0.0324£'i.ec- The correlated pairs are prepared by first cooling 
independent atoms in an external harmonic potential with the ground-state half-width of 
aE = 5a (frequency of 1 kHz ~ 30 nK). When the linear external potential is switched on, 
the atoms start moving in the direction of decreasing potential energy. Figure which 
captures the situation at three consecutive times, shows that unpaired atoms (off-diagonal 
peaks) are displaced by a much longer distance than the diatoms (diagonal peaks). 

The paired atoms remaining in the lattice are then in the state ~ exp[— (j — 
jo)'^(i'^ /{'2cro)]\xf^)\xf^) wherein positions and momenta are correlated with the uncertain- 
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FIG. 10: Simulation of the EPR state preparation in an optical lattice with 25 sites, at three 
consecutive times: (a) Initially (t = 0), the atoms are cooled down to the external harmonic 
potential ground state, whereas the LIDDI is off. (b) At f = 1.4 x lO"'^ s LIDDI and the repulsive 
linear potential (with the slope 0.04 Ej-ec per lattice site) are on, whereas the harmonic potential 
is off. The diatoms are moving through the lattice very slowly in comparison to the single atoms, 
(c) At t = 2.16 X 10^^ s single atoms are ejected out of the lattice and the diatoms are separated 
out. 

ties Ax+ ~ aE/V^ and Ap^ ^ h/Ax+, respectively. At higher temperatures the atoms are 
not cooled to the ground state of the external potential and the momentum anti-correlation 
has the spread 



Ap+ ^ h/{V2aEtii-nh[h^/{2almkBT)]}. 
The parameter s of Eq. (jHT|) can then be estimated as 



(35) 



V2i 



tanh 





/ a ^ 


\2 E 
\ -'-'Tec 






' kBT_ 



(36) 



This relation enables us to select the optimal external harmonic potential (specified here by 
(T^;) such that the parameter s is maximized, for a given temperature T. 

The small effective mass of unpaired atoms allows us to cool them to temperatures higher 
than that corresponding to the bottom of the diatomic band. The price is, however, that 
most of the atoms are discarded and only a small fraction of ~ a/cTE will remain in the 
bound diatom state. The different behavior of the paired and unpaired atoms in a periodic 
potential is a sparse-lattice analogy of the transition from Mott-insulator to a superfluid 
state in the fully occupied lattice, recently observed in Ref. 
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The two-particle joint position distribution of the ground state is a chain of peaks of half- 
width a separated by a that are located along the line X2 = xi (Fig. [TT)). The corresponding 
joint momentum distribution spreads over an area of half- width h/{2a) and consists of ridges 
in the direction p2 = —pi, that are separated by 2'Kh/a, and have the half- width TTh/{Na) 
for a lattice of sites (Fig. IT^ . 



VI. MEASUREMENTS 



After preparing the system in the EPR state, how can one can test its properties ex- 
perimentally? To this end we may increase the lattice potential Uq, switch off the field 
inducing the LIDDI, and separate the two lattices by changing the laser-beam angles. By 
increasing Uq, the atoms lose their hopping ability and their quantum state is "frozen" with 
a large effective mass: the bandwidth Vb decreases exponentially with Uq and the effective 
mass increases exponentially, so that the atoms become too "heavy" to move. One has then 
enough time to perform measurements on each atom: 

a) The atomic position can be measured by detecting its resonance fluorescence. After 
finding the site occupied by atom 1, one can infer the position of atom 2. If this inference is 



(a) 


o 








© 


X2/a 

< 












® 




o 





-6 -4 -2 2 4 6 
Xi/a 





2.5 


II 


2 




1.5 








1 








0.5 








(b) 



-1 




-0.5 





X2/a 



0.5 



FIG. 11: (o) Joint probability distribution of the positions of two lithium atoms in adjacent optical 
lattices, prepared in a diatom state as specified in the text, using the ground state of the external 
harmonic potential with half-width aE = 6a and temperature of 10 nK. (b) Position probability of 
atom 2 in the state above, conditional on atom 1 being measured at site (full line). Dashed line: 
Gaussian approximation of the Wannier function with the half- width a = 0.14a. 
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confirmed in a large ensemble of measurements, it would suggest that there is an "element 
of reality" Q corresponding to the position of particle 2. An example of the conditional 
probability of position of particle 2 after measuring the position of particle 1 is given in Fig. 

b) The momentum can be measured by switching off the x-lattice potential of the atom 
(thus bringing it back to its "normal" mass m). The distance traversed by the atom during 
a fixed time is proportional to its momentum. An example of the conditional probability of 
the momentum of particle 2 given the momentum of particle 1 is shown in Fig. ^1 

c) One can test the EPR correlations between the atomic ensembles occupying the two 
lattices, testing large number of pairs in a single run. The correlations in x and anti- 
correlations in p would be observed by matching the distribution histograms measured on 
atoms from the two lattices. 




FIG. 12: (a) Joint probability distribution of the atomic momenta in the aforementioned state 
with T = 100 nK. (6) Conditional probability of the momentum of atom 2, given that the mo- 
mentum of atom 1 has been measured (the measured value pi = piu is indicated by an arrow)for 
lithium diatoms prepared as in the text. The dashed line corresponds to the marginal probability 
distribution of momentum p2 irrespective of the momentum of atom 1 at the temperature T = 100 
nK. The half-width of each peak is equal to 1/s of Eq. 
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VII. CONCLUSIONS 



We have discussed a scheme which can be used to prepare a translationally entangled pair 
of massive particles in a state analogous to the original EPR state A novel element of 
the present scheme is the extension of the EPR correlations to account for lattice-diffraction 
effects. Their momentum and position correlations principally differ from those of free 
particles [Eqs. EPRstateDelta, EPRstateGauss]: due to the lattice periodicity, the position 
and momentum distributions have generally a multi-peak structure. 

The realization of the proposed scheme is expected to be based on the adaptation of 
existing techniques (optical trapping, cooling, controlled dipole-dipole interaction), to the 
requirements spelled out in Sec. |3 and IVll The most important ingredient of the scheme is 
the manipulation of the effective mass, for EPR-pairs preparation (by separating the "light" 
unpaired atoms from the "heavy" diatoms) and for their detection (by "freezing" the atoms 
in their initial state so that their EPR correlations are preserved long enough). 

n 

One may envision extensions of the present approach to matter teleportation and 
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Such extensions may 



quantum computation based on continuous variables |25l . 
involve the coupling of entangled atomic ensembles in optical lattices by photons carrying 
quantum information. 
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